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Abstract. Active galactic nuclei (AGN) are variable on a wide range of 
timescales, though relatively few systematic variability surveys have been con- 
ducted. Previous broad-band (both spectral and temporal) variability surveys of 
AGN are limited in their temporal and spectral bandwidth, despite their promise 
for probing the central engine and black hole mass. We outline optimal proper- 
ties for variability studies and provide a brief summary of three new variability 
surveys, two (BATSS and DASCH) about to begin and the third (EXIST) pos- 
sible within the next decade, which will open new windows on the physics and 
fundamental properties of AGN. 



1. Introduction 

A number of AGN studies have been undertaken for the variability of both 
optical (e.g. Peterson, these proceedings) and X-ray fluxes (e.g. Uttley, these 
proceedings) to enable constraints on the mass of the central supermassive black 
hole (SMBH) and to constrain emission models for the central engine, accretion 
disk and jet(s). Most of these have been with targeted narrow-field telescopes 
with observational cadence and total duration necessarily limited. The ideal 
AGN variability survey (AVS) would have the following properties, Pn: 

1. AVS-P1: broad sky coverage, fi, to maximize the number of AGN observed 
and to enable rare classes of variable objects, and low duty cycle events, 
to be found; 

2. AVS-P2: long total survey duration, D, of survey observations which each 
detect minimum source flux S m i n on timescale r Q and with fractional un- 
certainty in flux 6S To and enable variations on timescales from a maximum 
Tmax ~D/2 down to a minimum (median) r m i n ~ 0.7D/N to be measured 
from N randomly sampled observations; and 

3. AVS-P3: large total number of measurements, N, to enable measures of 
source variability on timescales r n ~ (n— l)r m ^ n /2, where n = 2 ...N, and 
fractional variability sensitivity improves as SS Tn ~ 8S To /n ^ 

Previous AGN variability studies have typically met only AVS-P3, and then 
usually with relatively short duration D. Broad-field (AVS-P1) AGN variability 
surveys are almost unknown, though the Swift/BAT survey (Markwardt et al 
2005) with ~70% sky coverage per day has begun to open up this domain. Here 
we outline three new surveys that will each extend one or more of these AVS 
properties. 
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2. New AGN variability surveys 

Time domain studies (e.g. PanStarrs and LSST) will unleash new constraints 
on AGN parameters and models by having one or more of the broad properties 
listed above. Already, optical and hard X-ray timing studies for full-sky AGN 
samples are beginning, and a far-reaching X-ray/7-ray temporal-spectral survey 
could EXIST as an Einstein Probe in NASA's Beyond Einstein Program. 

2.1. DASCH: Optical variability on scales At - lOd -(50-100y) 

Over the past 3y, we have developed an astronomical plate digitizer (Simcoe et 
al 2006) that is some 100X faster than any previously built in order to make 
available the Digital Access to a Sky Century from Harvard (DASCH). This 
will make possible (grant funding or a donor permitting...) the digitization and 
on-line access to the full images and derived photometry of Harvard's unique 
collection of some 600,000 astronomical images (all at least 5° x 7° ) of the full 
northern and southern sky from c. 1880 - 1985. Astrometry to <1" is derived 
from WCS solutions for each scan and photometry from SExtractor isophotal 
analysis fits calibrated locally from the >3000 GSC2.2 stars (B ~8-15) typically 
on each plate (Laycock et al 2007). An example light curve (Fig. 1) for a random 
star in the open cluster M44, used for development of photometric analysis 
software, shows the <,0.1mag (rms) photometry possible which for this dataset 
from 5 different plate series over 88y. Analysis of an initial sample of 15 bright (B 
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Figure 1. Partial lighturve (B mag vs. plate no.) for star (#399=GSC1395- 
2447) in the Open Cluster M44 measured on 245 plates (1890 - 1978). Even 
smaller rms errors are likely when annular calibration (on GSC2.2) is used vs. 
the full plate average done here since psf variations are then included. 

<,15) PG QSOs, starting with 3C273, is planned for a pilot study of variability 
power density spectra (PDS) to explore PDS break timescales as a measure 
of SMBH mass. With ^1000 plates for any given object randomly observed 
over ^100y, the median sampling time is ^25d and so the possible variability 
timescale range is r max /r m i n ~ 50y/25d ~ 730. Allowing for >,3 timescale 
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measures above a PDS break to determine r^, the corresponding SMBH mass 
range could be constrained over dynamic range of ^240. 



2.2. BATSS: Hard X-ray variability on scales r < lOOs-ld 

We have also initiated a "BAT Slew Survey", BATSS (Copete et al 2007) using 
the BAT hard X-ray imager (Barthelmy et al 2005) on Swift to analyze "event- 
mode" data from the ^60 slews (^l-2min each) that Swift performs each day 
to slew on/off pointed targets. Whereas BAT pointings cover some ^70% of the 
sky each day, adding in the slews increases sky coverage to nearly 100% as well 
as provides the only high time resolution data (apart from GRBs) since BAT 
pointing data is binned on 5-7min timescales. BATSS will thus provide AVS- 
P1,P3 and be particularly well suited to detect rare, bright AGN flares such 
as the extreme Blazar events from PKS2155-304 for which Swift /XRT /BAT 
coverage did not quite overlap with the 8 and 17Crab(!) TeV flares reported 
by HESS (Foschini et al 2007). Although the XRT spectra indicate that the 
synchrotron spectral break for this Blazar is below the BAT band, the BAT 
Transient Monitor (Krimm 2007) clearly does see flare variability from others - 

e.g. Blazar Mrk 421 (Fig. 2) for which extreme flares could be seen by BATSS. 
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Figure 2. Swift/BAT 2y lightcurve for Blazar Mrk421 showing ~100mCrab 
flares that might be resolvable (and still brighter) by BATSS. 



2.3. EXIST: Ultimate hard X-ray variability on scales r < 10s-5y 

The best prospects to optimize properties AVS-P1-P3 are with the EXIST mis- 
sion (htt p://exist.gsfc.nasa.gov/] ), proposed as the Black Hole Finder Probe in 
NASA's Beyond Einstein Program. EXIST images the full sky 3-600 keV each 
95min orbit with two large area and field of view (FoV) coded aperture tele- 
scopes (Grindlay 2005 and Grindlay et al 2007). With daily full-sky flux sen- 
sitivity S m i n ^lmCrab (comparable to Swift/BAT in ly) due to nearly 20% 



Figure 3. Simulated ly EXIST logN-logS for 21° x 21° field of view of one 
sub-telescope containing 331 AGN and 4° x 2° zoom image including brightest 
(4mCrab) source and 12 other AGN with fluxes from ~0.05 - 0.28mCrab. 

continuous coverage on every source enabled by continuous scanning with its 
large FoV and total area, EXIST would detect and study >3 x 10 4 AGN full 
sky. Each is located to <11" (90% confidence radius) which allows unambigu- 
ous host galaxy identification for its 0.05mCrab (=5 x 10 -13 cgs, 40-80 keV) 5a 
survey threshold sources. A simulated ly survey image and logN-logS is shown 
in Fig. 3 with normalization from Treister and Urry (2005) including obscured 
AGN. Thus, ^300 AGN (full sky) can be measured on timescales r m i n — Id or 
-1000 AGN with r min = 6d. For a 5y mission, AVS-P3 gives N = 1800 and 300 
timescales, respectively, to constrain the PDS and and thus SMBH mass. 

3. Prospects for SMBH mass, demographics and evolution 

Given the vs. mass scaling found by Uttley et al (2002) from RXTE data 
on several AGN, EXIST and possibly DASCH can constrain SMBH masses in 
AGN over a range -10 7 " 95 M , with >lO 9 M traced out to z <3. Non-AGN 
SMBHs (e.g. SgrA*) can be identified by their tidal disruption of main sequence 
stars and accompanying hard X-ray flares detected by EXIST and optical flares 
detected by DASCH for non-obscured systems. Together with extreme Blazar 
events from BATSS, the prospects for new AGN variability surveys are timely. 
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